Fig.lc : Growth of Power spectrum of density fluctuations. Spectrum and epochs 
are same as in Fig. la. 

Fig. Id : Power spectrum for potential is plotted against linearly extrapolated fluctu- 
ations in density. The dashed curve corresponds to the expected behaviour in linear 
theory. The dashed curve has been plotted for a = 2. Spectrum and epochs are same 



as in 


Fig. 


la. 






Fig. 


le : 


Power spect 


rum for gravitational force 


. All other details are same as Fig. 


Id. 










Fig. 


2a : 


Same as Fig 


. Id but for n = 0. 




Fig. 


2b : 


: Same as Fig 


. le but for n = 0. 




Fig. 


3a : 


Same as Fig 


. Id but for n = — 1. 




Fig. 


3b : 


: Same as Fig 


. le but for n = — 1. 




Fig. 


3c : 


: Nonlinear power spectrum for density 


is plotted against linearly extrapo- 



lated power spectrum at the same scale for n = — 1. Dashed line corresponds to linear 
behaviour. Notice the almost linear evolution of the nonlinear density contrast upto 
A(k) ~ 20. 

Fig. 4a : Same as Fig. Id but for n = —2. 
Fig. 4b : Same as Fig. le but for n = —2. 

Fig. 4c : Same as Fig. 3c but for n = —2. Here the nonlinear density contrast 
deviates from the linear density contrast in the quasilinear regime. However, in the 
nonlinear regime, both the curves are parallel to each other. 

Fig. 5a : The power spectrum of density perturbations is plotted on the n — A plane 
for CDM spectrum with T = 0.2. Curves have been plotted for five epochs, a = 0.1, 



0.25, 0.5, 1, 2. 












Fig. 


5b 


: Same 


as 


Fif 


5. 5a but for model with T = 


0.35 


Fig. 


5c : 


Same 


as 


Fig 


;. 5a but for model with T = 


0.5. 


Fig. 


5d 


: Same 


as 


Fif 


5. 5a but for model with T = 


0.7. 


Fig. 


6a : 


Same 


as 


Fig 


;.la but for CDM, T 


= 0.2. 




Fig. 


6b 


: Same 


as 


Fif 


5.1b but for CDM, T 


= 0.2. 




Fig. 


6c : 


Same 


as 


Fig 


;.la but for CDM, T 


= 0.35. 




Fig. 


6d 


: Same 


as 


Fif 


5.1b but for CDM, T 


= 0.35. 




Fig. 


6e : 


Same 


as 


Fig 


;.la but for CDM, T 


= 0.5. 




Fig. 


6f : 


Same 


as 


Fig 


.lb but for CDM, T : 


= 0.5. 




Fig. 


6g : 


Same 


as 


Fig 


;.la but for CDM, T 


= 0.7. 




Fig. 


6h 


: Same 


as 


Fif 


5.1b but for CDM, T 


= 0.7. 
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by plotting the power spectrum on the n — A plane [ See fig. 5]. It is seen that for small T, 
the spectrum is flat and the physically relevant scales enter the nonlinear regime almost 
simulateneously, killing the conspiracy of indices and amplitude to some extent. However, 
for large T, power spectrum is steeper and the region with n = —2 is well into the nonlinear 
regime when scales with n = — 1 are in the quasilinear regime. On the other hand, spectra 
with large T have stronger nonlinearities as compared to models with small T, leading to 
smaller variations in ^{k). Another feature worth noting in Fig. 5 is that for these spectra, 
a large range of scales develop a similar index, indicating power law correlations in that 
range of scales. 

The results show that in hierarchical clustering models, two different factors combine 
together in deciding the evolution of the gravitational potential. Smaller scales reach the 
quasilinear and nonlinear phase at an earlier epoch and if their indices are appropriate, 
their growth rate becomes similar to that of linear evolution. On the other hand the 
curvature of the spectrum helps to bring in the right scale at the right amplitude. In real 
life, of course, all these are of only approximate validity; but it illustrates a very interesting 
dynamical aspect of the nonlinear gravitational clustering. By and large, the gravitational 
potential seems to be a rugged entity in the hierarchical models changing less than what 
might have been anticipated by naive arguments. 

JSB would like to thank CSIR India for the Senior Research Fellowship. 
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Figure Captions 

Fig. la : Power spectrum for gravitational potential is plotted as a function of scale. 
Curves have been plotted for Ps{k) = Ak n with re = 1 at five epochs, a = 0.1, 0.25, 
0.5, 1,2. At the left end, lower curves correspond to later epochs. 
Fig. lb : Power spectrum for gravitational force. Spectrum and epochs are same as 
in Fig. la. 
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This shows that if n = — 1, then £ oc a 2 even in the quasilinear regime. It is not clear how 
significant this result is and whether it is possible to provide a simple analytic argument 
to prove the same. [This issue is currently under investigation.] 

Figures 4a and 4b give the results for n = —2. We see that there are deviations from 
linear theory in quasilinear regime but nonlinear evolution is similar to the linear one. 
This, of course, is a consequence of the scaling mentioned above. This result, taken in 
isolation, is not too interesting because: (i) the evolution in the quasilinear regime changes 
the amplitude, thereby making concrete predictions difficult and (ii) the result could be 
obtained trivially from the stable clustering argument applied to self-similar evolution. 
Figure 4c shows analogue of figure 3c for this spectrum. 

But the above two results [for n = — 1 and —2] taken together lead to an interesting 
"conspiracy" for a certain class of power spectra. To see this consider a power spectrum 
which has an index —2 in the nonlinear regime and —1 in the qusilinear regime. In the next 
instance, all the scales will grow in proportion to the expansion factor all the way from 
the strongly nonlinear regime to the linear regime. If we now smooth the spectrum and 
arrange its curvature such that the above condition is also [at least approximately] satisfied 
in the next instance, then the linear evolution will last longer. For such a spectrum, the 
gravitational potential will evolve very little. 

Of course, this requires conspiracy between the slope and amplitude of the power 
spectrum at different scales. Surprisingly enough, CDM like spectra do have such a con- 
spiracy built in to some extent. Figure 5c shows the local index of the COBE normalised 
CDM spectrum as a function of the density contrast, as the spectrum evolves. In order 
to maintain approximate constancy of the gravitational potential, the index should stay 
around —2 for A 2 J> 200 and around —1 for 1 <S A 2 <S 200. Figures 5a, 5b and 5d give the 
corresponding curves for different values of T for a CDM like spectra parametrised as 



Here 6, c and d are functions of T. [We are following the parametrisation of Efstathiou et 
al., 1992]. It is usually claimed that T = 0.2 is a good fit to observations. We see that the 
conspiracy is stronger for smaller values of T as the nonlinearities reached for these spectra 
are much lower than for spectra with high T. The evolution of C and D for these spectra is 
shown in figures 6a, 6b, 6c and 6d. In each figure we give the curves for a = .1, .25, .5, 1, 2. 
All the spectra are normalised to COBE. We see that the potential can be treated to be 
approximately constant [say, to 10 percent accuracy] up to scales of 2h~ 1 Mpc at a = 1.0 
in the standard CDM. For the "best-fit" case of T = 0.2 there is no appreciable change in 
A^(k) upto a = 1. However, if we compare models at the same level of nonlinearity at a 
scale of lh~ 1 Mpc then spectra with higher value of Y fare better. This can be understood 



P(k) 



Ak 



(16) 



[(1 + bk + {ckfl 2 + (dk) 2 Y] 



2jv 
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The self similar scaling of these curves can be understood in the following manner. 
Consider A^(k) } it is related to the power spectrum of density fluctuations as 



= ^r- (12 » 

We know that for power law spectra, power spectrum of density fluctuations must be a 
universal function of [k/k n [(a)] where k n i(a) is the scale that is going nonlinear at the 
epoch a ; say, the scale at which A(k) is unity. This is also true of the linearly evolved 
power spectrum. If we can express A^(k) in terms of a self similar function of [k / k n i(a)] 
and a, we would have obtained the relevant scaling. By rewriting (12) we obtain 

A^(k) = -Jj^A(k/k nl )(^j cxa- 2 ^f(k/k nl ). (13) 

We have used the fact that k n \ oc a n+3 in deriving the scaling. Similar relation can also 
be derived for A g (k) } 

A g (k) = J^A(k/k nl )(^j cxa- 2 ^f(k/k nl ). (14) 

Here we notice an interesting fact : for n = 1, A^(k) is a nonevolving function of [k / k n [(a)]. 
Similarly for n = — 1, A g (k) is a time independent function of [k / k n i(a)]. 

Figures 2a and 2b give curves corresponding to figure Id and le for n = spectrum; 
since this contains the relevant information we have not plotted graphs corresponding 
to figures la, lb and lc in this case. The two scales, corresponding to quasilinear and 
nonlinear regime, are apparent in the graph. 

We get our first surprise in the case of n = — 1 spectrum, shown in figures 3a and 3b. 
We now see that the evolution follows linear result even in quasilinear regime, all the way 
up to A(k) ~ 20! In figure 3b, A g (k) has been plotted for all the epochs mentioned above 
but as A g (k) is a time independent function of A^(A;), the curves fall on top of each other. 
The corresponding evolution of density contrast is shown in figure 3c by plotting A(a, k) 
vs. Ai(a, k) at the same scale. It is seen that the nonlinear evolution follows the linear 
evolution closely upto A(a, k) ~ 20. [As far as the authors know, this fact has not been 
specifically noted previously in the literature.] This fact shows that n = — 1 spectrum 
[corresponding to the isothermal density profile of p oc x~ 2 ] is rather special; for this 
spectrum, linear theory results have a validity which is beyond its legitimate domain. 

It is, of course, possible to demonstrate this more explicitly using the fitting law in 
(7). For a powerlaw spectrum with index n, we can explicitly calculate the scaling relation 
with a using this and we find that 

i/ \ ( -, // , ,\ r.n + 1 + 1)0+3) — 

Z(a,x) _ J (,7)i/("+4) a -^, 1 + JU (forl<e< 200 ) n 

a x K ' { (11.7y /(n+b) a »+sx n+a (for 200 < £ ) 
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concentrate on the scalar 

D(a,x) = (g(a,y + x).g(a,y)) (10) 

The power spectrum associated with D(a,x) is related to P(a } k) by a k 2 factor. It is 
easy to see that D (a, x) and £ (a, x) are related by the equation 

1 dD (a, x) 1 £ (a, x) 

This can be integrated readily to obtain D(a } x). A comparison with the linear correlation 
Dl(o,,x) is useful for estimating validity of quasilinear approximations at a given scale 
and epoch. Largest scale for which the change is above tolerance limit provides a natural 
criterion for validity of approximation schemes like frozen potential. Such a criterion can 
also be used to select a smoothing scale for truncated Zeldovich approximation. The 
comments made above regarding C are also applicable to D. 

3. Conclusions 

We shall now present the results of our analysis. We first consider pure power law spectra 
for density with linear power spectrum P (k) oc k n and n = 1,0,-1,-2. For sake of 
uniformity all of these are normalised to give erg = 1 at a = 1 with a gaussian window 
function. Initial power spectra for density and potential are Ps (a, k) = Aa 2 k n ; P^ (a, k) = 
Ak n ~ A with 

4 7 r 2 (8// i )( 3 +") 
" r[(3 + n)/2] 

Figures la and lb give the evolution of A^(k) and A g (k) for n = 1 spectrum. Curves 
are plotted for five epochs a = 0.1, 0.25, 0.5, 1, 2. With our normalisation we expect the 
scale of / = 8h~ 1 Mpc to go nonlinear around a = 1. We see that at the nonlinear end 
A^(k) and A g (k) change by a substantial amount, though the change here is less than the 
corresponding change in density correlations. Figure lc gives the evolution of A(k) during 
the same epochs. 

The results for power law spectra, of course, can be presented in a much more mean- 
ingful way. Our basic aim is to look at changes in potential in the nonlinear regime. For 
this, we should plot A^, and A g as a function of A^(A;), the linearly evolved power spec- 
trum for density fluctuations. Since there is no intrinsic scale in this problem, we expect 
this curve to evolve in a self-similar manner. Figures Id and le show these plots for n = 1 
spectrum. It can be shown below that A^(k) is a time independent function of A^k) (see 
below), this implies that curves corresponding to different epochs will coincide, as shown 
in figure Id. 
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2. The Gravitational Potential and Force 

The results in (6) and (7) relate the mean correlation funtion for density perturbations 
in the linear and nonlinear regimes. To obtain the information about the gravitational 
potential from this result, the most direct approach for studying the evolution of gravi- 
tational potential is the following: The correlation function C (a, x) for the dimensionless 
gravitational potential ?/>(a,x) = (2/3) Hq~ 2 Qq~ 1 Lp (a, x) defined as 



C(a,x) = (ip(a,y + x) ?/> (a, y)) y (8) 
is related to the matter correlation function £(a, x) by the relation 



V 2 V 2 C 



1 d ( 2 d 



1 d ( 9 dC (a, x) 



(9) 



where we have assumed that C (a, x) = C (a,x). The above equation follows immediately 
from the fact that the power spectra for density and potential are related by a factor 
k 4 in the fourier space. Given a correlation function in linear theory one can use (6) or 
(7) to obtain the nonlinear correlation function; integrating equation (9) we can find the 
evolution of correlations in gravitational potential. While integrating this equation, one 
can use the linear theory results at very large scales to obtain a unique solution. By 
fourier transforming the correlation function one can obtain the power spectrum of the 
gravitational potential. 

It is, however, possible to simplify the above procedure by using the following fact : For 
spectra which vary sufficiently smoothly, the quantity A 2 (a,k) = [k 3 P(k) (a, k) /27r 2 ] at 
k ~ x~ x behaves in a manner very similar to £ (a, x). So, instead of finding the correlation 
function of potential by integrating (9) and then obtaining the power spectra by fourier 
transforming, one can directly use the ansatz (6) in the fourier space to obtain the power 
spectrum of gravitational potential. [A similar relation has been suggested by Peacock 
and Dodds [Peacock and Dodds, 1994] for evolution of A 2 . As this relation is almost 
identical to (6) for = 1, we use (6) to obtain the nonlinear power spectrum. In their 
paper Peacock and Dodds make cautionary remarks about the applicability of (6) for 
spectra with n < — 1. However the key results described in this paper depend only on the 
characteristic behaviour in the quasilinear and nonlinear regime as described in (7), which 
are seen in simulations of all spectra of physical interest.] 

Such a simplification, of course, needs to be justified. We have used both methods 
and ascertained that the deviations are not significant in regimes of interest. 

In addition to the gravitational potential we also investigate the evolution of the 
gradient of the potential, g(a,x) = — which is more directly relevant in determining 
the dynamics. The correlation function for g(a,x) is a six-component object but we shall 
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where x = /). In other words, one can find the nonlinear evolution of the correlation 

function and related quantities from the linear theory using the above ansatz. It is, how- 
ever, possible to produce a still simpler fitting function which captures the essence of this 
formula: [Bagla and Padmanabhan, 1993] 



As we shall see, this simple form allows one to understand clearly the effects which are 
operating at different scales. It is apparent that evolution of density contrast is generically 
characterised by two different scales: At about £1, ~ 1 the deviations from linear theory 
begin to manifest itself. The second scale occurs around £1, ~ 6 [with £ 200] when we 
expect virialised structures to form and separate out from the overall dynamics. These 
two scales are characterised by the different slopes in (7). For the sake of definiteness we 
shall call the three regimes linear, quasilinear and nonlinear. 

Before proceeding further, we would like to make some comments regarding the ac- 
curacy of the above ansatz. This relation was first obtained by Hamilton et al., based 
on N-body simulation data of Efsthathiou et al [ Hamilton et al., 1991; Efstathiou et al., 
1988]. In the original paper, it was claimed that the relationship is reasonably good and 
valid at all epochs and for all spectra. Recently, this relation has been tested by more 
accurate simulations with larger dynamic range [Padmanabhan et al., 1995]. The results 
of this paper suggests that the "universality" is only approximately valid. The form and 
asymptotic limit of h has a weak spectrum dependence and the asymptotic value is lower 
for spectra with more small scale power. However, this analysis also shows that the devia- 
tions from the universality are small and are at about 20% level. To this level of accuracy, 
one can use the above ansatz. The accuracy can, of course, be easily improved by fitting a 
more accurate [spectrum dependent] curve to the results. We shall not do so, since the key 
idea of this paper is only to develop a simple physical picture of growth of the gravitational 
potential. 




(for < 1-2, C< 1-2) 

(for 1.2 < < 6.5, 1.2 < £ < 195) 

(for 6.5 < £ L , 195 < £) 



(7) 
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this paper. 

In principle, the question of evolution of gravitational potential can be settled in 
a straightforward manner by running suitable N-body simulations. However, such an 
approach does not provide one with an intutive understanding of the results. Because 
of this reason we shall follow a more indirect route in this paper which, it turns out, 
produces results that are equivalent to the N-body simulations (to the accuracy we want). 
This approach is based on the observation that numerical simulations suggest a simple 
relationship between the mean relative pair velocities of particles v (a, x) and the mean 
correlation function £ (a, x) at the same epoch and scale [Hamilton et al., 1991]. The mean 
correlation function £ (a, x) is defined as 

X 

£(a,x) = ^ j ' i{a,y)y 2 dy (2) 
o 

Let h(a } x) = [— v(a,x) /ax] be the dimensionless relative pair velocity at the comoving 
scale x at the epoch a. Then, following Nityananda and Padmanabhan [Nityananda and 
Padmanabhan, 1994] we postulate that h (a, x) is a universal function of £ (a, x) evaluated 
at the same epoch and scale; that is 



h (a, x) = Q [£ (a, x)] 
where the function Q can be well approximated by 



(3) 



HO 



2g 
(3+?) 

560 

exp 



-0.23 



40 



(for f < 38) 

(for 38 < £ < 100 ) 

(for 100 < 



(4) 



It can be shown that [Nityananda and Padmanabhan, 1994] this result allows one to express 
the true mean correlation function £ (a, x) in terms of the correlation function in the linear 
theory £l ( a ? as 



(2 rZ( a > x ) 
3 



dq 



h(q)(l + q) 



(5) 



where / 



[i + e>, 



1/3 



For the functional form in (4) the relation between £ (a, x) 



and £l ( a 7 can be well approximated by [Hamilton et al., 1991] 



C(a,x) 



x + 0.358x 3 + 0.0236x 6 
1 + 0.0134x 3 + 0.00202x 9 / 2 



(6) 
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1. Introduction 

The driving force behind the formation of large scale structures in the universe is the 
gravitational field produced by inhomogeneities. Overdense regions accrete matter at the 
expense of underdense regions allowing inhomogeneities in the universe to grow. At scales 
with L <C H -1 where H = (a /a) it is essentially this process which allows small initial 
inhomogeneities to grow into structures like galaxies, clusters etc. Further, observations 
suggest that the material content of the universe is dominated by dark matter which 
may be considered to be made of collisionless elementary particles. In that case, the 
gravitational force is dominated by these particles and, to first approximation, we can 
ignore complications arising from baryonic physics. The evolution of inhomogeneities is 
then governed purely by the gravitational force. 

It is, therefore, intersting to ask the question: How does the gravitaional potential in 
the universe evolve as structures form ? In addition to the inherent academic interest, this 
question has two other facets to it. 

To begin with, we note that the gravitational potential <p (a, x) due to inhomogeneities 
is governed by the Poisson equation 



where 8 (a, x) is the density contrast and pi, is the background density. [ This equation is 
valid in the limit of L <C H -1 where L is the scale we are interested in. Such a description 
is quite adequate for the study of nonlinear structure formation.] When the perturbations 
are small and the linear theory is applicable, 8 oc a if fio = 1 and the universe is matter 
dominated. [If Ho ^ 1, the gravitational potential undergoes slow evolution during the 
linear regime that can be computed using linear theory. In this paper we shall concentrate 
on the fio = 1 model.] Hence the nontrivial evolution of gravitational potential takes place 
only in the nonlinear phase. Because of this reason the study of growth of gravitational 
potential offers a direct diagnostic of the nonlinear dynamics. 

The second aspect is the following: Sometime back the authors have suggested an 
approximation scheme for the study of gravitational dynamics called "frozen potential 
approximation" [Bagla and Padmanabhan, 1994; this was also suggested independently by 
Brainerd et al, 1993]. Comparison of this approximation scheme with N-body results shows 
that it works well in the quasilinear regime and provides a reasonably accurate information 
about velocities. In fact, for certain class of spectra this scheme works unexpectedly well, 
showing that the average effect of growth of gravitational potential is not as important as 
one would have naively imagined. This requires an explanation. In particular, one would 
like to know whether there exists a class of spectra for which gravitational potential does 
not evolve significantly even in the nonlinear regime. We shall address these questions in 
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Abstract 

We study the evolution of the power spectrum of gravitational potential during the nonlinear 
clustering in an ft = 1 matter dominated phase. N-body simulations suggest that the potential 
does not evolve in time even in the quasilinear phase for an n = — 1 power spectrum. For 
n = —2, the potential evolves in the quasilinear phase but not in the extreme nonlinear phase. 
Becauase of these facts, the evolution of the gravitational potential in spectra like CDM is less 
than what would have been expected naively. We discuss a class of CDM like models in which 
such an interesting conspiracy between the amplitude and local index occurs. 
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